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Abstract 

In this work a result of existence and uniqueness for a plane cavity driven steady 
flow is deduced using an analytical method for the resolution of a linear partial 
differential problem on a triangular domain. The solution admits a symbolic 
expression based on integration over the domain. Some examples of flow are 
computed and graphed. In particular, it is shown a realistic example of a shear- 
forced flow with two eddies, usually computed only by numerical methods. The 
mathematical techniques used for the demonstration of the main result are ele- 
mentary. 
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1 A theorem of existence and uniqueness 

Let fl be (the interior of) a triangular domain in R^, the cartesian {x, y} plane, 
with vertices O — (0,0), A — (2a, 0), B = (a, a), where a is a positive real 
number. Note that the triangle OAB is rectangular and OB = BA. Let be 
f{x,y) G C'^{fl,R). We want to resolve the differential problem 

- dL^ + dyy<f> = f in^^, = on on (1) 

for a function (j){x,y), 4> S C^(0,]R). Note that the partial differential equa- 
tion —9^2.0 -I- dyycj) = f admits a general solution of the form (see [5] or [lOj) 
4>{x, y) = g{—x + y) + h{x + y) + <j)o{x, y), where g and h are arbitrary real func- 
tions and (j)o is a particular solution of the equation. But this general expression 
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is not very useful for applying the boundary condition cf) — on d^. It is more 
interesting and instructive the following direct method. 
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Figure 1: The domain in the xy -plane. 

Consider the differential operator —d^^ + dyy written as {dx+dy){—dx+dy), 
and consider a linear transformation rule for cartesian coordinates X = ax + by, 
Y — cx + dy. If we want to have 2dx — dx +dy and 2dY = —dx + dy, using the 
chain rule it must hea = b= l = d= l and c — ~1, that is 

X ^x + y, Y = -x + y (2) 
The transformation is invertible: 

2x ^ X -Y, 2y ^ X + Y (3) 

With the notation ^{X,Y) — <p{x{X,Y),y{X,Y)) and analogous for /, the 
differential equation —31x4' + 9yy4' — f becomes 

4dj,YHX,Y) ^ F{X,Y) (4) 

Note that the transformation ([2]) is a 45°-rotation and a v^-dilation of the 
plane {x,y}. Also, the boundary condition doesn't change: $ = on dil (for 
simplicity we use for the domain in the plane {X, Y} the same symbol used 
for the plane {x,y}). For example, ^{X,—X) = 0(a;,O) = 0. Therefore, the 
differential problem ([1]) becomes 

AdxY^ ^ F infl, $ = on 9f7 (5) 

We remark the fact that the operator 49^y is the canonical form of the dif- 
ferential operator — + dyy, which has the lines y = x and y = —x as 
characteristic curves ([S] or [lOj). 

Now we want to discuss the resolution of the differential problem (O . 
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Figure 2: The domain in the XY -plane. 
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Figure 3: Domain of integration. 



Let P = {X,Y) be a point in the interior of the domain fl. Then we can 
construct the polygon E using segments parallel to X and Y axes (see Figl3]). 
Note that M = {-Y,Y), R = (-^,0), Q = (2a, 0), S = (2a, -X), N = 
{X,—X). Using the identity 29^y — {dxdy + dydx): from the differential 
equation it follows that 
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2 [ [dxdY<^{X,Y) + dYdx^{X,Y)]dXdY ^ [ F{X,Y)dXdY (6) 
Now apply the Green theorem ([5] or [TT]) to the first integral: 

/ {dxdY<i> + dydx^jdXdY = [ {dy^ dY - dx^ dX) (7) 

It is now simple to calculate the line integral along the edges of the polygon E 
(note that the boundary must be walked in counter-clockwise sense): 

/ (ay$ dY - dx<^> dX) = 
= -2$(P) + 2$(iV) - 2$(S') + 2$(Q) - 2$(i?) + 2$(M) (8) 
So we have 

- 2$(P) + 2^{N) - 2$(S') + 2$(Q) - 2$(i?) + 2$(M) = 

= / FdXdY (9) 

Now apply the boundary condition $|ao = 0: it follows that 

$(P) = $(X,r) = -i /" F{t,s)dtds (10) 

with the consequence that, if the point P{X, Y) lies on the boundary of 51, that 
is if P = 7\f = N or P = S or P = R, the function P must satisfy the necessary 
condition 

0=1 F{t, s) dt ds V(X, Y) e dn (11) 
"'s(jf,y) 

It is easy to see that previous condition can be written in a more explicit fashion: 

/ / F{t, s) dsdt = yX e [0, 2a] (12) 
Jx J-x 

Therefore we have shown that a solution to differential problem ([5]), and hence 
to ([1]), exists if and only if F satisfies condition p2| . Also, formula (fTO|) is an 
analytical expression for a solution. Note that, denoted by T the point {X,0), 
the integral can be divided into the two integrals defined on the two simple 
rectangles PTRM and NSQT. 

Now we discuss uniqueness of solution. Suppose to have two solutions <&! 
and $2 for the problem Then $ = $i -$2 is a function such that dxy^ — 
y{X, Y) E n and ^\gn = 0. Note that we can write 
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dy [dx^^ = 2 dx^ aiy$ = (13) 
Applying the Green theorem to domain S for the expression dy [dx'^^ , we have 

[dx'^fdX = 

IdT, 

= [ [dx^fdX+ [ [dx<i>fdX+ [ [dx^fdX (14) 

Jn Jq Jm 

Using integration by parts, the following identity holds: 

j [dx'^f = $ 9x4" - j d^x^ dX (15) 
Hence, being $|go = 0, the second integral in is null, therefore 

/ [dx<^fdX+ / [dx<i>fdX = Q (16) 
JN Jm 

The two integrals are evaluated in the same sense of the integration path, so 
that = along the segments NS and MP, therefore $(P) = $(A/). But 
$(M) = 0, being M e dil, so for a generic point P — {X, Y) we have $(-P) = 0. 
The solutions <i>i and $2 are identical. 

We have shown the following result (remember that X = x + y, Y ^ —x + y): 

Theorem 1 Let f be a real function of C'^{n,M.) such that 

^t-s t + s\ 
f(^,-^]dsdt = 0VXe[0,2a] (17) 



X J-x 
Then the differential problem 

- dl^<p + d^y(l> ^ f inn, (f) = ondn (18) 

has one and only one solution in the space C^(r2,R). The solution is given by 
the formula 

H^, y)= - I f f f ^) ds dt - 

^ .2a M ft-S t + s\ 

-9 / fi^^^jd^dt (19) 

^ J x+y J—x—y \ ^ ^ / 
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2 An application: cavity driven flows 



In this section we discuss an application of previous theorem to a problem of 
two-dimensional cavity driven flow, that is a plane flow confined in a cavity and 
induced by the stress due to a primary flow external to the cavity (see [3]). 
This phenomenon has great importance in scientific research (see e.g. [7]) and 
technological applications. Assume that the cavity has the shape of the triangle 
OAB of FiglT] in the xy-plane. Stress due to the primary flow acts on the 
horizontal edge OA. We suppose that the fluid is ncwtonian and incompressible, 
that is plane stress T and plane strain-rate D tensors are linked by the formula 
(see m 

T = 2^]D) (20) 

where /i is the dynamic viscosity and 23ij ~ {djVi + diVj) (see f6|), where 
'v={vi,V2) = {vxTVy) is the flow velocity field. Plane incompressible flows admit 
a stream function ( 6 ), that is a function ^'(x, such that 

u = Vx= v^vy^ -dx'i' (21) 

Therefore, a plane newtonian incompressible flow is described by the partial 
differential equation 

^dl,^ + dyy^^-T,y (22) 

M 

In the next of the paper we suppose to know the analytical expression of T^y 
and we try to find a solution of for a stream function "if such that 5'|ao = 0. 
This boundary condition is usual for plane incompressible flow (see [B] and [2]), 
but in the case of a cavity driven flow it (or an analogous 5'|go = const) has an 
important physical meaning. In fact, if 5'|ao = 0, then dfl is a level curve for 
VP, therefore at each point of the boundary V^* is orthogonal to the tangent of 
the boundary itself ([H]). But V4' = {—v,u), which is orthogonal to the flow 
velocity field (u, v). Therefore, at each point of dil, the geometrical tangent and 
the velocity field are parallel, that is the flow is conflned into the cavity fl. 
Applying theorem ([T]), it can be stated that if Txy € C''(ri,M), then there is 
a unique stream function e C'^{fl,W) solving the linear equation (p2l) with 
boundary condition 'i'\on — 0. Note that 1 consider a nonlinear problem about 
cavity driven flow where uniqueness can fail. 

We consider at first the more simple analytical form for a possible stress: 

T^y = /i(ciy -I- C2) (23) 

In this case, along the horizontal edge OA {y = 0) of the cavity the stress 
is constant. From theorem ([1]), a solution to our differential problem exists if 
the function ciy + C2 satisfies the condition pT|) . It is easy to show that the 
condition is satisfied for all X £ [0, 2a] if and only if 2c2 = — aci. Note that in 
this case the stress has expression 
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T,y = AiC2— y+1 (24) 



a 

and for ?/ = I it changes its sign. So flow can recirculate. We make the choice 
C2 = —8a, so that ci = 16. The differential problem to solve is 

- dl^(j3 + dly<f) = IQy - 8a in 17, = on dn (25) 

From formula ^9]) . using the transformation rule X = x + y,Y = —x + y the 
solution to ([25]) can be computed by 

^{x, y) = - 4: / {t + s - a) ds dt - 

- 4 / / (t + s-a) dsdt (26) 

J x+y J —x — y 

which gives the expression 

4'(x, y) = 2y^ - 2x^y - 4aj/^ + 4aa;y (27) 

for the stream function of the flow. The velocity field is {dy'$, —dx'^) = (— 2a;^ + 
6y^+Aax—8ay, Axy—Aay). It is interesting to find the points where the velocity is 
null. Solving the algebraic system (— 2a;^ + 6t/^ + 4aa:: — Say, 4a;?;— 4aj/) = (0,0), 
we find as expected the three vertices (0,0), (2a, 0) and (a, a), and also the 
interior point (§, a) which is the center of the recirculation gyre (see Figl4]). 




Figure 4: Flow path-lines in the case a=l. 



Now we consider a more interesting case. Let the stress be described by a 
sinusoidal expression of the form 

Txy = Afi cos(fcy) (28) 

with A and k real numbers. Using ^T7\ . it is easy to show that if we suppose 
q = 0, then 



7 



k — m— , 
a 



TO = 2?! + 1, rt e N 



(29) 



is the condition for existence and uniqueness of a flow in the triangular cavity. 
Consider m = 1. By integration p9p . the analytic form of the stream function, 
solution of the differential problem, is 



2Aa^ 
97r2 



cos ( j + cos - y) 



2cos" ( J(a^ + y) 



(30) 



and FigO shows some path lines, where one primary central eddy and three sec- 
ondary eddies are present. It is also interesting to draw the graph of u = dy'^ 
for X = a and y variable in the range [0, a]: there are two values of y, not equal 
to a, for which m = 0. One of the two values is equal to the value where u = 
in the previous case of the linear stress (see Figl6]) . 




Figure 5: Flow path-lines in the case of the sinusoidal stress, with a=l, A=5. 




Figure 6: Comparison of u = dy^i in the linear (thin line) and sinusoidal (thick 
line) case. 



A more realistic example is based on a stress with analytical expression of 
the form 



4 

T,, -/i (31) 

m,n— 

Applying condition (I17|) for the computation of the coefficients a„i_„, a possible 
stream function is 



^{x, y) = (2y3 _ 2x^y - 4ay^ + 4axy){y - lOOx^ - a) i^y + -x - -aj (32) 

The horizontal component u of the velocity field, along the x-axes, is a 5-order 
a:-polynomial whose graph is shown on FiglT] The stress acts on the horizontal 
segment of the triangular domain as a variable shear of positive sign. 

u(x,0) 
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Figure 7: Profile ofu{x,0) — dy'i'{x,0) in the case of a variable shear stress. 
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Figure 8: Flow path-lines in the case of a variable shear stress, with a— I. 

The resulting flow path- lines (see FigE]) show the presence of a primary gyre, and 
of a secondary gyre near the vertex opposite to the edge subjected to external 
stress. This image is similar to a corrisponding picture (fig. 2(b)) in [2], where 
stream function is computed by numerical method. 
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